Free magnetized knots of parity-violating deconfined matter in heavy-ion collisions 
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We show that the local parity violation in the quark-gluon plasma supports existence of free 
(meta)stable knots of deconfined hot quark matter stabilized by superstrong magnetic fields. The 
magnetic field in the knots resembles the spheromak plasma state of the magnetic confinement 
approach to nuclear fusion. The size of the knot is quantized, being inversely proportional to 
the chiral conductivity of the quark-gluon plasma. The parity symmetry is broken inside the knot. 
Particles produced in the decays of the knots have unusual azimuthal distribution and specific flavor 
content. We argue that these knots may be created in noncentral heavy-ion collisions. 

PACS numbers: 25.75.-q,25.75.Nq,11.30.Er 



MOTIVATION 

Heavy-ion collisions may generate very strong mag- 
netic fields with the magnitude of the order of the QCD 
scale. The lowest bound of the maximal strength of the 
magnetic field B that can emerge in noncentral collisions 
of Pb-Pb ions at the Large Hadron Collider (LHC) at 
CERN was estimated in [1| to be 
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where e = |e| is the absolute value of the electron charge 
and m,r * 140 MeV is the pion mass. This superstrong 
long-ranged magnetic field acts as an external field on 
electrically charged quarks and antiquaries that consti- 
tute the hot and dense quark-gluon matter created in 
the collision. 

Theoretically, the magnetic field of the strength (JTJ 
may lead to various nonperturbative effects such that (i) 
a modification of the phase diagram of QCD via inequiv- 
alent shifts of both the chiral [2j and deconfining finite- 
temperature transitions (ii) an enhancement of the 
chiral condensate Q via an enforcement of the dynamical 
chiral symmetry breaking H; (iii) an anomaly- mediated 
creation of a stack of parallel domain walls made of neu- 
tral pions in a dense plasma etc. 

Strikingly, the superstrong hadron-scale magnetic field 
has already exposed itself experimentally. Recently, the 
STAR collaboration has found possible signatures of 
parity violation in locally macroscopic domains of hot, 
strongly interacting matter created in heavy-ion colli- 
sions at the RHIC experimental facility Q ■ Such domains 
may arise due to a nonvanishing expectation value of 
the gluon topological charge that breaks both parity (P) 
and charge-parity (CP) symmetries. In the hot quark- 
gluon plasma the spontaneous generation of the parity- 
violating domains - suggested theoretically in Ref. Q - 
may be revealed with the help of a very strong magnetic 
field which is created in a noncentral heavy-ion collision 



as well. The key point is that in the parity-violating 
background the external magnetic field leads to genera- 
tion of an electric current. This unusual phenomenon - 
used in the experiment [7j to reveal the local parity vi- 
olations via electric charge fluctuations - is now known 
as the Chiral Magnetic Effect (CME) [l,[loj]. There exist 
also a numerical evidence of existence of the CME based 
on the lattice simulations of QCD [ll| . 

The essence of the CME may be formulated in the 
following way: if a magnetic field B is applied to a system 
characterized by an asymmetry between the number of 
right- and left-handed fermions (quarks), then an electric 
current j is induced along the magnetic field axis [lfj| : 
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a x B. 



(2) 



The external electric field is assumed to be absent, E = 0. 

In ([2]) the quantity a x is the chiral conductivity (loj 
and 
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is the conserved electromagnetic four-current of the 
quarks. The vector j represents the spatial components 
of ([3]). In Eq. ^ e q and tpf are the electric charge and 
the field of the / th quark flavor, respectively, and 7 M is 
the Dirac matrix. The left-right quark asymmetry in the 
quark-gluon plasma - that is needed for a realization of 
the CME ([2|) - may be induced by various parity-odd 
(topological) effects [lj|. 

The parity-odd law @ can be contrasted with the ordi- 
nary Ohm's law, j — oE, which is certainly a parity-even 
feature of the system. Here a is the conventional conduc- 
tivity, and the external magnetic field is assumed to be 
absent, B = 0. 

We show that the presence of the chiral conductiv- 
ity (2} - induced by topological QCD effects - may lead 
to existence of (meta)stable knots made of the magnetic 
field that is frozen in the hot quark-gluon plasma in 
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a magnetohydrodynamic equilibrium. In general, knot- 
ted configurations are of wide interest in many areas of 
physics. One can encounter them in particle physics fl2j ]. 
fluid dynamics 13 1 , plasma physics [14j , condensed mat- 
ter physics , and classical field theory [Jjj], E3 ■ Knot- 
ted light solutions to the vacuum Maxwell equations are 
known theoretically (lij , and very recently they were dis- 
covered experimentall y in laser beams in a complemen- 
tary approach of Ref. [19j . Here we show that free knot- 
ted objets can be created in heavy-ion collisions as well. 



KNOT SOLUTIONS IN PARITY-ODD PLASMA 

The classical Maxwell equations of motion are 

d^=j v , d^ v = 0, (4) 

where F^ v is the electromagnetic field strength tensor, 
= (l/2)e^ al3 F a f i is its dual, and j M is the electric 
current density ([3]). It is convenient to rewrite Eqs. (|4]) 
in terms of electric E and magnetic B fields: 



VxB 


8E -. 
~ ~dt ~ J '' 


(5) 




V-E = p, 
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V-B = Q, 


(8) 



where p = jo is the electric charge density. 

We are interested in static electrically neutral solu- 
tions, and therefore we set 



dB 

E = 0, p = 0, - 



0. 



(9) 



The quark-gluon plasma is conducting, and therefore 
the Maxwell equations should be supplemented with the 
laws, that account for effects of the medium on the prop- 
agation of the electric current. In the absence of the bulk 
electric field © this transport property is determined by 
the CME, Eq. ©. 

The requirement (|9]) is automatically consistent with 
Eqs. ((6|) and ((7|), while Eqs. §5§ and dHJ) give us a simpler 
system of equations: 



VxB = a x B . 
V ■ B = 0. 



(10) 
(11) 



We are looking for solutions of Eqs. (TTUl) and (fTTT) in- 
side a ball of a certain radius R. As we discuss below, the 
quark matter inside the ball, r < R should be in a decon- 
fined phase while the state outside the ball corresponds 
to the confinement phase. Therefore, the electric current 
j in a stable configuration must always be tangent to the 
boundary of the ball, j ■ r = 0. Here the unit vector in 
the r direction, f, represents the outward-normal vector 



to the boundary of the ball. Then Eq. ([2]) gives us the 
following boundary condition: 



B(r) 



= 0. 



(12) 



The system of equations (fTU|) and (lll| is very well 
known because it plays an important role in usual plasma 
physics. The corresponding solutions are usually called 
Chandrasekhar-Kendall states following the authors who 
stressed importance of these solutions for astrophysical 
plasmas [20| . In the theory of nuclear fusion plasmas, 
the eigenvectors of (fT0|) and (jlip are known as Taylor 



states 
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Such states are ultimate minimum-energy 
configurations in the plasma equilibrium. 

Mathematically, Eqs. ([10]), CO} and ((T2j) define an 
eigensystem of the curl operator "Vx". Usually, an 
eigenvalue problem is formulated in a fixed domain of 
space, so that the unknowns are eigenvalues. In our case 
the eigenvalue in Eq. (|10p is fixed by the known chi- 



ral conductivity a x , which is assumed, for the sake of 
simplicity, to be coordinate-independent. The unknown 
quantity is the radius of the ball, which enters the system 
of equations via the boundary condition (|12l) . 

In the spherical domain the curl operator has a discrete 
set of quantized eigenvalues and eig envectors, that were 
described in details in Refs. [20j, |2lj, |23j . Adopting these 
results to our notations, we find that Eqs. (fTTT)) . (fTT|) . (fT2|) 
have solutions in spherical balls of the quantized radii, 
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where /cj? is the fcth positive zero of the Bessel function 
J„ + i/2(k). In our paper we are interested in the knot 
solution with the smallest possible radius 
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k = = 4.49341 



(14) 



In the spherical coordinates the corresponding magnetic 
field configuration is 

B{r, 9, <p) = B [u{a x r, 6)r + v(<r x r, 9) 9 

+w(a x r,9)<p], (15) 

for r ^ R , and B(r, 9,ip)=0 for r > R . In Eq. ^B) 
2 
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u(f , 0) = 73 (sin £ - £ cos £) cos ( 



(16) 
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e ) = 73 ( sil ^ - e cos£ - £ 2 sin£) sin0 , (17) 



w{£,6) = |2 (sin£ -£ cos £) sinS, 



(18) 



and 9 and tp are the unit vectors in the azimuthal and 
polar directions, respectively. The azimuthal angle, 9 6 
(0, 7r) is inclination from the z axis, while the polar angle, 
ip G (— 7r,7r), operates in the xy plane. The coordinate- 
independent factor B in Eq. (IT51) is a free parameter. 
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The structure of the magnetic field of the solution (|15j) - 
(IT8| is rather complicated. The solution consist of an 
infinite number of nested (deformed) tori, and the de- 
gree of the deformation varies with the torus size. The 
magnetic field lines wrap around each torus and pro- 
duce multiple self-linkages, what is a typical feature of 
all Chandrasekhar-Kendall states 
the solution (HSl) - llHl) as "knot" . 



18j. Therefore we call 



as "knot" . The structure of typ- 
ical magnetic lines and electric currents in the knot is 
illustrated in Fig. [1] 





FIG. 1. Typical structures formed by the lines of the mag- 
netic field and, equivalently, by the electric currents. The 
lines wind around different nested tori inside the same knot 
solution (|15|l - (|18|) . The figures are described in the text. 

It is convenient to label each structure by specify- 
ing the maximal size r max of the (deformed) torus in 
units of the overall knot size Rq, Eq. (fT4|) . In the 
largest structures the field wind around two tori simulta- 
neously [Fig.QJa), r max /i? = 0.96]. As the size becomes 
smaller, the structure of the magnetic fields gets a con- 
ventional torus shape [Fig. fjjb), r max /i?o = 0.91]. Yet 
another torus shape can be found at lower radii [Fig.QJc), 
fmnx/Ro — 0.77]. At particular radii the magnetic field 
lines form well-known knotted structures. For exam- 
ple, at r max /i?o = 0.88566, the magnetic lines are trefoil 
knots, which wrap around a torus, Fig. Hid). 

Notice, that despite the knot occupies an ideally spher- 
ical region, its interior is only axially (not spherically) 
symmetric. Moreover, the knot is a parity-odd object: 
its reflected image cannot be superimposed on itself. A 
mirror image of a left-handed knot is a right-handed knot. 

In an astrophysical applications the knot (| 1 5 [) - (|T5|) was 
used to model magnetic field patterns of the Crab Neb- 
ula [13] ■ At much higher energies, in the nuclear fu- 
sion concept, the solution (fl~5 |) - ([T8l) is considered as a 



relaxed minimum energy state of a plasma in a sphero- 
mak device (Le. , inside a spherical shell that confines the 
plasma) [HI, [2^. In general, the solutions of the system of 
equations (fTOf , (jTTJ) and (fT2|) are the "natural end config- 
urations" in a magnetically dominated plasma, because 
the magnetic turbulence drives such plasma towards the 
spheromak state 0, HB} . We suggest that similar stabi- 
lized configurations may arise in the quark-gluon plasma 
created in heavy-ion collisions. 



BASIC PHYSICAL FEATURES OF KNOTS 

Phase. The interior of the knot should be in a decon- 
fined state with unbroken chiral symmetry. The decon- 
finement is needed for the quarks to flow freely inside 
the knot, so that they should be able to produce stabi- 
lizing magnetic fields. The chirally restored symmetry 
is also very important because in this state the chiral 
condensate is zero. A presence of the chiral condensate 
would wash out the imbalance between the left-handed 
and the right-handed quarks, leading to vanishing chiral 
conductivity [lOj], a x = 0. The very existence of the knot 
is guaranteed by a nonvanishing value of a x , and there- 
fore the knot cannot exist if its interior is in the chirally 
broken phase. The chirally restored and simultaneously 
deconfined phase is expected to be produced at the RHIC 
facility and at the ALICE/LHC experiment [26] . so that 
the knots have a chance to be created there. 

Generally, one may expect that the cold exterior out- 
side the knot should provide an additional factor that 
stabilizes the knot. The nonzero tension of the interface 
between the phases would squeeze the knot increasing 
the energy density of the magnetic field, and generat- 
ing a backreaction force from the inside. A qualitative 
illustration of the knot structure is presented in Fig. [2] 




FIG. 2. A free hot knotted solution ([T5j)-(IT8j) in the cold QCD 
vacuum. The tori that host the magnetic lines of Figs. Dla)- 
(c), and some of the magnetic lines are shown. 



The electric charge of the knot is zero. Other impor- 
tant parameters of the knot may be expressed via two 
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free parameters of the solution (fT5j) - (TT8l) : the magnetic 
field scale Bo and the chiral conductivity a x . 

The knot radius Rq is given by Eq. (|14[) . 

The total magnetic energy of the knot is 



£ mag n = \f d 3 rB 2 (r) = C E ^L = 5 2 i? 3 , ( I.!) ! 



Bn Cj. 

■x ,? o 



where Ce = (47t/3)ko sin 2 no = 17.9337 

The total magnetic moment of the knot is 
1 - -/ ~ B a Cm 



M 



df-rxj = C A/ — I 



3 B i?3 . 5 , (20) 



where C M = ?r 2 [l - (1 + k§/2) cos/c ] = 33.6592 . . . and 
z is the unit vector along the z-axis. 

The flavor content of the knot is constrained by the re- 
quirement of stability of the knot with respect to direct 
quark-antiquark annihilations. This implies that the in- 
terior of the knot may consist of either u, d and s quarks, 
or u, d and 5 antiquarks. A mixed interior, say, u, d and 
s is less favorable energetically. The presence of heav- 
ier quarks - with the masses heavier than a few T c - is 
dynamically suppressed. 

The condition of the local electric neutrality [here rif = 
n f (x) is the density of quarks of / th flavor] , 



0, 



2e 



sets the constraint on the total quark numbers Nf. 

2N u = N d + N s with N f = I d 3 rn f (r) . (22) 



A similar relation holds for the antiquarks in a knot 
made of antimatter. Due to relatively large mass of 
the strange quark we expect that the strange fraction, 
p s = N S /(N U + N d + N s ), should be low, /3 S < 1. Thus, 
the flavor composition of the knot's interior should be 
similar to the one of a neutron star, N d — 2N U . 



The minimal baryon charge of the knot, 



Qi 



lJ2 N f = l( N " + N i- 



N, 



(23) 



can be estimated as follows. The local electric current is 

(24) 



f 



2e e e 

-7^n u v u - -n d v d - -n s v s , 



where Vf is the speed of the quarks of / th flavor. A mini- 
mal value of the baryonic charge is reached in the case of 
the maximal "charge-parity polarization" , i.e. when all 



positively charged quarks (u) are, say, left-handed while 
the negatively-charged quarks (d and s) are all right- 
handed. In this case the CME-driven ([2j current takes its 
maximum value at a fixed baryon charge. In the strong 
external magnetic field positively charged u quarks and 
negatively charged d and s quarks move in opposite direc- 
tions. Ignoring a contribution of the strange quarks (i.e., 
setting P s =0), one gets v u — —v d for the light quarks 
which travel almost with the speed of light, \v u ,d\ = 1- 
Using (f22|) . we get the absolute value of the current: 

j(r) = |/(r)| - ^n u (r) = ^n d (r) . (25) 

Equation ([2]) then allows us to express the quark's densi- 
ties n Uid via the magnetic field of the solution (TT5|) -([T^ | . 
Then the baryon charge of the knot (f2"3"]) is: 
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Bq _ Cq b r, 



eat 



(26) 

where Cq — 80.4854 . . .. The knots made of antimatter 
have a negative baryonic charge, Qb = —\Qb\- 

In a real physical environment of heavy-ion collisions 
the approximate characteristics of the knot can be es- 
timated as follows. At high temperature T the chiral 
conductivity is (loj : 



3g 2 



2tt 2 T 2 + ^ 2 /tt 2 ' 



— 1 

Q = — 
N f 



2 - — y, 
\ j, 



f 



(27) 



where Nf is the total number of participating quark fla- 
vors, q 2 is the average squared electric charge of these fla- 
vors, fi is the quark chemical potential, and 715 = c$- fm~ 3 
is the chiral density given by the difference between the 
number of left- and right-handed quarks per unit volume. 
Here C5 = 2NfQ w , where Q w is the topological winding 
number (accumulated) per fm 3 . 

For a qualitative estimation of the knot's parameters 
in Pb-Pb collisions at ALICE, we take Nf = 2 (so that 
q 2 = 5e 2 /18) and we assume that scale of the finite 
chemical potential is the same as the temperature scale, 
h/tt ~ T. Thus, slightly above the phase transition, 
T - 200 MeV, one gets a x - 0.4 c 5 MeV. Then the knot 
radius (H} is Rq ~ (2 • 10 3 /c 5 )fm. Typically @, [13, 
C5 = (1 — 10) so that the typical radius of the knot is 
as small as R ~ 200 fm. This length is large compared 
to a geometrical scale (~10fm) of a typical heavy-ion 
collision, so that a creation of such large knot is unlikely. 
However, very large fluctuations of the topological charge 
in particular events are not excluded. Thus, in our qual- 
itative estimation we set Q w ~ 50 per fm 3 , that gives us 
C5 ~ 200 and a x = 80 MeV. As for the scale of the mag- 
netic field Bo, we take a fraction of the peak value (fT|), 
e£>o ~ m 2 . Then we get the radius (fT4"|) . the total mag- 
netic energy (fl"9|). the total magnetic moment (f20| in 



5 



terms of the nuclear magneton /Ujv = e/(2m p ), and the 
minimal baryon charge (|23[) of the knot, respectively: 



|M| 



10 fm, 

2.5- IOVjv, 



niagii 

Ob 



150 GeV, 

2.5 -10 3 . (28) 



The radius of the knot Rq is of the order of the size of 
the lead nucleus. The magnetic energy of the knot -Emagn 
is also a natural quantity being a bit smaller than the 
rest energy of the Pb nucleus. The corresponding energy 



density, e e 



-Emagn/(47ri?g/3) ~ 35MeV/fm 3 is much 



smaller than the non-equilibrium energy density regime 
of ALICE, £ ~ (1 - 1000) MeV/fm 3 [26]. The magnetic 
moment of the knot is rather large in nuclear units, while 
\M\ ~ 15 fiB hi terms of the Bohr magneton \ib- 

The baryon charge Qb of the knot (|28l) is six time 
larger than the total number (~400) of nuclei in the 
colliding lead ions. However, the total multiplicity of 
produced particles in Pb-Pb collisions is estimated to be 
in the range of thousands 26j, so that the large baryon 
charge Qb of the knot may be within experimental reach. 

The lifetime of the knot is difficult to estimate from 
general considerations since a proper analysis should in- 
volve thermal magnetodynamic calculations in a so-far- 
unexplored parity-odd backgrounds. We hope that this 
analysis can be done using a numerical approach. 



GEOMETRY OF THE KNOT DECAY 

The knot may be observed via the particles produced 
by its decay. Besides the specific flavor content (the u and 
d quarks should appear in proportions 1:2), the flow of 
produced particles is rather unusual as we will see below. 
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FIG. 3. The contour plot of the strength of the electric cur- 
rents q — \j(r)\/(o'xBo) in the transverse xz plane (deter- 
mined by y — 0) of the knot. The contour lines (from outer 
to inner) correspond to g = 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.65. 



Although the knot occupies the spherical region, the 
knot's interior has a lower (axial) symmetry. The spatial 
distribution of the magnetic energy density has a form 
similar to the one of an oblate (disk-shaped) spheroid. 
The short axis of the spheroid is parallel to the magnetic 
moment (|20[) of the knot (z-axis in our notations), while 
the short axes are located in the plane of the knot (xy- 
plane). A contour plot of the strength of the electric 
currents in a transverse xz plane of the knot is shown in 
Fig. G3 The standard eccentricity of the electric currents 
in the knot is negative: 



(z 2 - x 2 ) 3 

{Z 2 +X 2 )j 



-0.168. 



(29) 



where (. . .)j = J dxdz . . . j(x, z). In the contrary, an 
initial geometry of a noncentral heavy-ion collision has 
an almond shape characterized by a positive eccentricity. 

A decay of a knot should produce particles that domi- 
nantly follow the current distribution of the quark's cur- 
rent in the equilibrium state of the knot. Using the ex- 
plicit knot solution ()15|) -(fT8 |) one can calculate the dis- 
tribution of the decay products 



d 2 N 

dcf) dip 



d 3 r|j(r)|<5(n(0,^)-j(^)) , (30) 



where j (r ) = j (f)/ \j (r) | is the unit vector in the direction 
of the quark current at point r, and h{4>, <p) is the unit 
vector pointing in the direction determined by the polar 
angle (p and the new azimuthal angle (f> = ir/2 — 9. The 
angle (j) gives inclination from the xy plane (we follow 
the standard notations used in physics of the heavy-ion 
collisions). 

An explicit calculation reveals that the distribu- 
tion (I3U1) is independent of the polar angle <p, so that the 
decay of the knot is axially symmetric with respect to 
rotations around the z-axis. However, the azimuthal dis- 
tribution of the particles emitted from the decaying knot 
is quite nontrivial (Fig. |4j as it differs significantly from 
a typical azimuthal distribution of particles produced in 
heavy-ion collisions. 
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FIG. 4. Normalized azimuthal distribution of particles cre- 
ated in the knot decay (the solid blue line), and its elliptic 
component (the dashed green line). 
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The azimuthal particle distribution can be expanded 
into a Fourier series, 



27r d-ZV 

-tt = 1 + 2 K cos(fc^) + ™k sin(/c0)] , (31) 

^ k=l 



were JVo is the total number of the particles produced 
in the knot decay. Symmetries of the knot imply that 
v 2 k-i = vo 2 k = with k g N. 

The first five nonzero v- and ^-coefficients of the 
Fourier expansion (|31j) are summarized in Table |U The 



TABLE I. The Fourier flow coefficients of the particles pro- 
duced in the knot decay. 



n 1 2 


3 


4 


5 


v 2 n 0.2378 -0.1138 
w 2 n-i 0.001695 0.01619 


0.04425 
0.005640 


-0.02389 
-0.002061 


0.01505 
0.001731 



direct flow of the knot decay is zero, v\ — 0, while the 
elliptic flow is rather large, v 2 ~ 1/4. The elliptic har- 
monic, 1 + 2v 2 cos(20), is shown in Fig. 0] by the dashed 
green line. Notice, that the higher-order w-harmonics are 
not small. For example, the ratio of fourth and (squared) 
second cosine coefficients is negative and large in its ab- 
solute value, 



Vi/vl = -2.01 . 



(32) 



contrary to the case of the heavy-ion collisions where one 
expects Vi/v 2 « 0.5 
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In the contrast with the heavy-ion collisions, the sine- 
terms of the expansion f|3 1 1) are nonzero. This fact indi- 
cates that the knot decay is asymmetric with respect to 
the z — > — z flips: the distribution of the emitted parti- 
cles, say, above the xy plane does not coincide with the 
distribution that is observed below the plane. The z-axis 
asymmetry is consistent with the fact that the magnetic 
moment of the knot is nonzero. 

The form of the azimuthal distribution, Fig. |4l signals 
that the largest flux of the decay products is emitted a bit 
above the knot's plane. The maximal flow is concentrated 
around the cone with the azimuthal angle 



0.343 . 



20°. 



(33) 



The flow of the decay particles in the direction of the 
knot axis (cf> = 0,7r) is zero. 

Finally, we notice that it is likely that the plane of the 
knot may tend to coincide with the reaction plane of a 
noncentral heavy-ion collision that creates the knot. In- 
deed, the strongest magnetic field - created in the center 
of the collision - is perpendicular to the reaction plane 
while the strongest magnetic field of the knot (observed 
in its center, Fig. [3]) is perpendicular to the knot's plane, 
Fig. [5] Since the knot is formed via a reconnection pro- 
cess of the lines of the magnetic field created in the col- 
lision, the orientations of the magnetic field lines of the 



collision and of the knot are likely to be correlated. This 
fact, in turn, implies a correlation of the corresponding 
planes. 




x/R 



FIG. 5. The projection of the magnetic field lines on the 
y — cross-section of the knot. The lines have, in general, a 
y component (not shown). 



SUMMARY 

In the summary, we outline our vision of the life of a 
typical magnetized knot. 

Creation. A noncentral collision of two heavy ions cre- 
ates a dense hot medium (plasma) of quarks and gluons 
pierced by lines of a strong magnetic field. The mag- 
netic field - that is induced by moving charged ions and 
products of their collision - is perpendicular to the reac- 
tion plane. The dense quark-gluon matter expands and 
cools down, while (some of) the magnetic lines reconnect. 
Simultaneously, fluctuations of the topological charge in 
the plasma lead to emergence of the local parity- violating 
domains characterized by the imbalance between the left- 
and right-handed quarks. The violation of the parity 
symmetry in the background of the magnetic field leads 
to the chiral magnetic effect, and the electric current start 
to flow along the lines of the magnetic field. 

Stabilization. If the reconnected magnetic lines hap- 
pen to be knotted, then the electric current creates a 
self-supporting state, since the current flowing at the one 
segment of the knot supports the magnetic field and, con- 
sequently, the electric current at other segments. Mag- 
netic turbulence drives the plasma towards a minimum 
energy configuration that resembles the Taylor sphero- 
mak state in the magnetically confined plasmas of nu- 
clear fusion devices. Both the closeness of the knotted 
magnetic lines, the conservation of the electric current, 
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and the squeezing effect of the knot boundary contribute 
to the stability of the knot. The magnetic lines hold the 
hot quark plasma confined inside the knot. The size of 
the knot is expected to be of the order of some tens of 
fermi. Generally, the knot is an electrically neutral sys- 
tem with a neutron-star-like flavor content: one u quark 
comes with two d quarks with a possible small admixture 
of s quarks. 

Decay. As the knot cools down, the emerging chiral 
condensate starts to wash out the imbalance between the 
left- and right-handed fermions and eventually leads to 
the decay of the knot. We cannot estimate the lifetime of 
the knot using general arguments, but we expect that it 
should be much larger than a typical lifetime of a quark- 
gluon fireball created in the collision. Primary products 
of the knot decay are u and d quarks coming out with 
the specific ratio 1:2. The emitted particles have the 
unusual azimuthal distribution characterized by a signif- 
icant elliptic flow V2 ~ 1/4, relatively large (and some- 
times negative) higher harmonics (i.e., v^jv\ ~ —2), and 
axially-symmetric broad cone-shaped distribution of the 
emitted particles. The aperture of the emission cone is 
180° — 2</> max w 140°. The specific flavor content and az- 
imuthal distribution of the decay products may serve as 
good experimental signatures of the knots at RHIC and 
ALICE/LHC experiments. 

The author is grateful to P. Forgacs and A.J. Niemi 
for useful discussions. 



* On leave of absence from ITEP, Moscow, Russia. 
[1] V. Skokov, A. Y. Illarion ov and V. Toneev , Int. J. Mod. 

Phys. A 24, 5925 (2009) [arXiv:0907.1396] . 
[2] E. S. Fraga and A. J. Mi zher, Phys. Rev.D 78, 025016 

(2008) [arXiv:0804~1452] ; Nucl. Phys. A 820, 103C 

(2009) [arXiv:0810.3693] . 

[3] N. O. Agasian and S. M. Fedorov, Phys. Lett. B 663, 445 
(2008) [arXiv:0803.3l"56] . 

[4] I. A. Shushpano v and A. V. Smilga, Phys. Lett. B 
402, 351 (1997) jarXiv:hep-ph/9703201] ; T. D. Cohen, 
D. A. McGady and E. S. Werbos, Phys. Rev. C 76, 



055201 (2007) arXiv:0706.3208;. 
[5] V. P. Gusynin, V. A. Mir ansky and LA. Shovkov y, Nucl. 

Phys. B 462, 249 (1996) [arXiv:hep-ph/9509320| . 
[6] D. T. S on and M. A. Step hanov, Phys. Rev. D 77, 014021 

(2008) |arXiv:0710.1084] . 

[7] B. I. Abelev et al. [S TAR Collaboration ] , Phys. Rev. Lett. 

103, 251601 (2009) |arXiv:0909.1739] ; ImXiv:0909. 17171 

see also a Viewpoint of B. Miiller, Physics 2, 104 (2009). 
[8] D. Kharzeev, R. D. Pisarski and M. H. G. T ytgat, 

Phys. Rev. Lett. 81, 512 (1998) |arXiv:hep-ph/980422lj ; 

D. Kha rzeev and R. D. Pisarski , Phys. Rev. D 61, 111901 

(2000) |arXiv:hep-ph/9906401| . 
[9] D.E. Kharzeev, L.D. McLerran, and H.J. Warringa, Nucl. 

Phys. A 803, 227 (2008) |arXiv: 071 1.0950] , 
[10] K. Fukushima, D.E. Kharzeev, and H.J. Warringa, Phys. 

Rev. D 78, 074033 (2008) |arXiv:0808.3382] . 
[11] P. V. Buividovich et a l, Phys. Rev. D 80, 054503 

(2009) iarXiv:0907.0494| : M. Abramczyk, T. Blum, 
G. Petropoulos and R. Zhou. larXiv:0911.1348l [hep-lat], 

[12] L. D. Faddeev and A . J. Niemi, Nature 387, 58 (1997) 

[arXivihep th/9610193| . 
[13] H. K. Moff'att, J. Fluid Mech. 35, 117 (1969). 
[14] M. A. Berger, Plasma Phys. Control. Fusion 41, B167 

(1999); L. D. Fa ddeev and A. J. Niemi, P hys. Rev. Lett. 

85, 3416 (2000) [arXiv: physics /0003083] . 
[15] E. Babaev, Phys. Rey . Lett. 88, 177002 (2002) 

[arXiv:cond-mat/0106360] . 
[16] E. Radu and M. S. Volkov, Phys. Rept. 468, 101 (2008) 

[arXiv:0804.1357] . 
[17] R. A. Batty e and P. M. Sutcliffe, Phys. Rev. Lett. 81, 

4798 (1998) [arXiv:hep-th/9808129| . 
[18] A. F. Ranada, Lett. Math. Phys. 18, 97(1989); W. T. M. 

Irvine, D. Bouwmeester, Nature Phys. 4, 716 (2008). 
[19] M. R. Dennis, R. P. King, B. Jack, K. O'Holleran, 

M. J. Padgett, Nature Phys., 6, 118 (2010). 
[20] S. Chandrasekhar and P.C. Kendall, Astrophys. J. 126, 

457 (1957). 

[21] M. N. Rosenbluth and M. N. Bussak, Nucl. Fusion 19, 
489 (1979). 

[22] J.B. Taylor, Phys. Rev. Lett. 33, 1139 (1974). 

[23] J. Cantarella, D. DeTurck, H. Gluck, M. Teyteld, Phys. 

Plasmas 7, 2766 (2000). 
[24] L. Woltjer, Proc. Natl. Acad. Sci. USA 44, 489 (1958). 
[25] J.B. Taylor, Rev. Mod. Phys. 58, 741 (1986). 
[26] ALICE Collaboration, J. Phys. G: Nucl. Part. Phys. 30, 

1517 (2004). 

[27] C. Gombeaud and J. Y. Ollitr ault, Phys. Rev. C 81, 
014901 (2010) |arXiv:0907.4664] . 



